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Abstract 

Cosmologically motivated theories that explain small acceleration rate of the Uni- 
verse via modification of gravity at very large, horizon or super-horizon distances, 
can be tested by precision gravitational measurements at much shorter scales, such 
as the Earth-Moon distance. Contrary to the naive expectation the predicted correc- 
tions to the Einsteinian metric near gravitating sources are so significant that might 
fall within sensitivity of the proposed Lunar Ranging experiments. The key reason 
for such corrections is the van Dam-Veltman-Zakharov discontinuity present in lin- 
earized versions of all such theories, and its subsequent absence at the non-linear 
level ala Vainshtein. 



1 Generalities 



Recent observations suggest that the Universe is accelerating on the scales of the 
present cosmological horizon This indicates that, either there is a small vacuum 
energy (or some "effective" vacuum energy), or that conventional laws of Einstein 
gravity get modified at very large distances, imitating a small cosmological constant 
0,0,1^. The first possibility is unnatural in the view of quantum field theory^, 
since the required value of vacuum energy ~ lO^^^GeV^ is unstable under quantum 
corrections. This unnaturalness goes under the name of the Cosmological Constant 
Problem. 

In this respect the second approach, of modifying gravity in far infrared can be 
more promising since it is perturbatively stable under quantum corrections. The 
unnaturally small value of the vacuum energy is replaced by the idea that laws of 
conventional gravity break down at very large distances, beyond a certain crossover 
scale Tc- The value of is perturbatively-stable ie. it does not suffer from cut-off 
sensitive corrections experienced by vacuum energy. This modification at the level 
of linearized effective four-dimensional equation for the freely propagating metric 
fluctuations /i^jy(x) about the fiat background can be described in the following 
general way i,0,i|Ti: 

+ 4/(r.^V^) ) V = (1) 

c 

Here = V^V^ denotes four-dimensional d'Alambertian, and f{'r'1V^)/r1 is an 
operator that dominates over only for momenta g -C I/tc. One way to think 
of this modification is that the effective gravitational coupling (effective Newton's 
"constant" ) 
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becomes dependent on the wave-length. We will assume that the function f{r'^V^) 
is expandable in power series of V^. 

Depending on the precise structure, such models can be divided in the following 
three categories. The first category is the model 0, in which Newtonian gravity 
turns into the five dimensional l/r^-potential at distances r ^ rc. This effect is 
due to existence of an infinite-volume fiat extra dimension to which gravity "leaks" 
from the 3-brane, where the conventional particles live. So a brane-observer sees 
effectively four-dimensional theory of gravity, in which the effective 4D graviton is 
unstable, but with an arbitrarily large lifetime r ~ over which it decays into 
five-dimensional continuum of states. 

Although extra dimensions play crucial role in formulating the manifestly generally- 
covariant theory, for our purposes we can equally well use purely 4D language. For 
the brane observer interested in metric on the brane created by brane sources, there 
is an effective four-dimensional description in terms of a single four- dimensional 



^An alternative explanation could be provided by Anthropic approach 
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graviton hfj_^. Not surprisingly, this effective 4D theory is non-local, but non-local 
terms only dominate in far infrared. Moreover, if one is interested in estimates 
of sub-leading corrections to Einstein metric near the gravitating sources, the 4D 
language is simpler. At the linearised level this 4D theory is given by (|l|) with 

/(r^V^) = r,v^ (3) 

In this theory the smallness of cosmological constant is not explained. However, 
if one postulates that vacuum energy is zero due to some other reason, the model 
explains the accelerated expansion of the Universe without any need of dark energy 
The reason, as discovered in is that in p| the effective Friedmann equation 
gets modified by additional powers of Hubble parameter H, that dominate for very 
low curvatures (late times). In the leading order the resulting equation can be 
parameterized as [|],|^ 

«'-^ = j7T- (4) 



At late times, this equation has a self-accelerating cosmological branch with H = ^. 
Observation fix the crossover scale to be Tc ~ lO^^cm. Confronted with minimal 
models where dark energy is a pure vacuum energy, the above model also contains 
just one parameter (rc versus the vacuum energy), but the difference is that the 
value of Tc is insensitive to quantum loops. 

In two other classes of models 0,^, one can address the cosmological constant 
problem . Models of are based on higher-dimensional generalization of p[ with 
N > 1 extra dimensions P]. As a result, the feature of 4D Newtonian gravity 
turning into the high- dimensional 1/r^^'^-potential for r ^ Tc is shared by iV > 1 
theories. However, having more than one extra dimensions proved to be crucial for 
cosmological constant. In these theories the large 4D vacuum energy (S) curves the 
four- dimensional space only mildly resulting in an accelerated expansion, with the 
rate inversely proportional to a power of ^. In four-dimensional language, gravity 
is modified in such a way that sources of wave-length ^ gravitate very weakly. 
Thus although vacuum energy is huge, it actually does not curve the space and can 
lead to an observable small acceleration. 

Finally there is a third class [§|,|T^ , which does not refer to any underlying high- 
dimensional theory, but is rather based on infrared modification of gravity directly 
in four-dimensions. In this models, Geff is postulated to act as a high-pass filter, so 
that as in [0 sources of wavelength ^ rc gravitate extremely weakly. 

In the above models, value of Tc may of may not be fixed. For instance, in the 
models of [0 it is actually fixed to be Vc ~ 10^*^ cm , due to the current lower limits 
on Newtonian 1/r gravity. 

The common feature of all the above theories, is that beyond the crossover scale 
Tc Newtonian 1/r gravity gets substantially modified. But, at short distances (near 
the gravitating sources), all the prediction of Einstein gravity are reproduced up 
to some small corrections. These corrections are the central interest of the present 
work. 
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Given the fact that crossover scale is one could naively conclude that devi- 
ations are suppressed by powers of ^. If this were the case, then for Vc ~ lO^^cm, 
the corrections would be completely negligible at any distance of potential precision 
measurements. 

Interestingly, the story is very different. Recall that near any gravitating object 
there is a second (much smaller scale) in the problem, the gravitational radius of the 




body rg. It was shown in ||T^ that leading corrections to the gravitational potential 
^ in the model of are strongly enhanced by inverse powers Vg and are given by 

(5) 

The aim of the present work is to show that in a large class of theories that modify 
gravity beyond some horizon or super-horizon size distance Tc, there are corrections 
of the type [|13| that penetrate at much shorter scales r,, <C Tc, and could be po- 
tentially measurable. These corrections can be detected in precision gravitational 
measurements in systems that are much smaller than Tc, e.g., such as the Earth- 
Moon system. In most interesting cases, such as the expression (^, they are on 
the border line of existing measurements, and thus can be detected if the existing 
sensitivities are improved. 

The key reason for such an unexpected behavior is that the graviton in the 
above-discussed theories has extra polarization that also couples to the conserved 
energy-momentum source, and mediates a scalar-type force. As a result, at the 
linearized level gravity is of that scalar-tensor type. Therefore whenever linearized 
approximation is valid the predictions of the above theory differ from those of Ein- 
stein gravity by a finite amount, no matter how small is the ratio r/r^. 

This effect was originally pointed out in the framework of linearized theory of 
massive gravity , and goes under the name of van Dam - Veltman - Zakharov dis- 
continuity (vDVZ). It lead vDVZ to the conclusion that massive theories of gravity 
are ruled out. Later the similar effect was observed in the linearized version of the 
generally-covariant model of [Q]. 

The crucial point [0, however, is that discontinuity is an artifact of the lin- 
earized approximation, and is cured by non-linear corrections. In reality, the solu- 
tion is continuous in the limit l/r,, — > 0. This was originally suggested by Vainshtein 
in |TD|0 and later confirmed by explicit fully-nonlinear analysis in generally covari- 
ant model of PI, both in cosmological solutic 
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such as cosmic strings [0 and Schwarzschild ||13|,|l6|. These studies uncover the 
same persistent pattern. Near the sources the solutions are arbitrarily close to those 
of Einstein gravity, and continuous for large Vc- However, the corrections set in at 



^ On AdS space vDVZD can be absent already at the linear level, if graviton mass is smaller 
than the AdS curvature |^,^,|2l[. The roots of this phenomenon can again be attributed to the 
effect of |lo| . Such situation probably is not experimentally interesting interesting, and won't be 
discussed here. 
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distances, (or time scales) much smaller than r^, which makes them potentially ob- 
servable even if is very large (e.g. horizon size) (the relevance of this corrections 
for the orbit of Jupiter was pointed out in |TB[). This makes us think that analo- 
gous behavior must take place in any theory with infrared modification of gravity 
of the form in (|I]). The key message of our analysis is that in all such theories of 
interest vDVZ discontinuity indeed disappears at the non-linear level in accordance 



to 113, [IT], [12], |13| , and standard predictions of Einstein gravity are recovered near 
the gravitating sources, where non-linearities are important. But as in [|13| for any 
given localized source, there exists a distance for which the non-linearities are 
unimportant and thus the "wrong predictions" take over. 

We give a qualitative prescription to estimate the leading order corrections to the 
Einstein metric using the form of the function /(t'cV^) in theories interest. Using 
this prescription one can derive observational constraints on such theories. 

The rest of the paper is organized as follows. In section 2, we briefly summarize 
our theoretical results and in section 3 we discuss their relevance for observations, 
through the anomalous perihelion precession of orbits of planets. The reader inter- 
ested in purely the phenomelogical applications of our results can dismiss the rest 
of the paper which is devoted to more theoretical considerations. In section 4 we 
discuss the role of extra graviton polarization in creating the observable effects at 
intermediate scales in theories with infrared-modified gravity. In section 5 we briefly 
discuss the applicability range of our results. 



2 Results 

Here we shall briefly summarize our results. We will be interested in a class of 
generally-covariant theories, containing exclusively spin-2 states, for which the lin- 
earized perturbations about the flat space satisfy: 

( + l^J{rlV') ) V = - ^ [t,. - f ^..T:| + . . . (6) 

where the ellipsis stands for some derivative-dependent terms that vanish in con- 
volution with conserved sources. The over- all coefficient on r.h.s. was absorbed in 
redefinition of Mpi. In all theories of interest /3 7^ 1, as opposed to the case of 
Einstein theory of massless graviton with two polarizations, in which j3 = 1. Thus 
the theory exhibits the analog of vDVZ discontinuity. 
Then our analysis suggests the following. 

As suggested by Vainshtein, vDVZ is cured at the non-linear level. That is, the 
solution are continuous in the limit ^ 00 near the sources. For any gravitating 
source, with gravitational radius -C there is an important intermediate scale 
in the problem 

< r,^ < Tc (7) 
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where r^, is determined from the following equation 




(8) 

The various distance scales work as follows. For r <^ r* the metric produced by 
the source is nearly Schwarzchild, with the leading correction to the metric that can 
be estimated as 

This expression, for the particular form of the /-function given by (0), correctly 
reproduces the result of [|l^ given in (|^), found by solving the full nonlinear high- 
dimensional equation of |0. 

As will be discussed below, the above corrections can play the crucial role for 
testing theories in question in precision gravitational measurements. 

In the interval r^, <^ r -C Vc, the linearized approximation is valid. So gravity is 
still 1/r, but is of scalar-tensor type. Note that the term "scalar-tensor" only refers 
to the similarity in the tensor structure, rather than to existence of an independent 
spin-0 state in the theory. Instead the "scalar" admixture comes from the extra 
polarization of spin-2 graviton. 

Finally, for the distances r ^ Tc, gravity becomes unconventional, with the 
potential given by 



where q is the three-momentum. 



3 Anomalous Perihelion Precession 

The considered class of theories predict slight modification of the gravitational po- 
tential of a massive body at observed distances according to (|^). This modifications 
can be observable in the experiments that are sensitive to anomalous perihelion 
precession of planets. Let e be the fractional change of the gravitational potential 

where \E' = —GM/r is the Newtonian potential. The anomalous perihelion preces- 
sion (the perihelion advance per orbit due to gravity modification) is 

50 = vrr(r2(r"ie)')', (12) 

where ' = d/dr. 

Let us apply this to the model of 0. In this theory, [|1^ 

e = -V2r-h-^/h^'\ (13) 
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The numerical coefficient deserves some clarification. The above coefficient was 



derived in |jT3| on Minkowski background. However, non-linearities created by cos- 
mological expansion can further correct the coefficient. One would expect these 
corrections to scale as powers of TcH , where H is the observed value of the Hub- 
ble parameter. On the accelerated branch [Q, as it's obvious from (^), i7 ~ 1/rc 
and thus, one would expect the corrections to be of order one. Recently, a very 
interesting fact was pointed out by Lue and Starkman '^^^ that the cosmological 



background only affects the sign of the coefficient. The sign depends on the partic- 
ular cosmological branch. It is negative for the standard cosmological branch, and 
positive for the self-accelerated one. We will restrict ourselves to order of magnitude 
estimate, but the sign will be very important if the effect is found, since according 



to [O] it could give information about the cosmological branch. 
We get 

50 = (37r/4)e. (14) 

Numerically, the gravitational radius of the Earth is Tg = 0.886cm, the Earth-Moon 
distance is r = 3.84 x lO^'^cm, the gravity modification parameter that gives the 
observed acceleration without dark energy Tc = 6 Gpc. We get the theoretical 
precession 

50 = 1.4x10"^^ (15) 

This is to be compared to the accuracy of the precession measurement by the lunar 
laser ranging. Today the accuracy is = 2.4 x 10~^^ and no anomalous precession 



is detected at this accuracy [|Tj . In the future a tenfold improvement of the accuracy 
is expected |r8| . 

As noted in |T^, anomalous Martian precession is also of interest for testing 
with cosmologically intersting values of Tc. For Tc = 6 Gpc we get 50Afars = 4x 10~^^. 
This is to be compared to the accuracy of cr^j^j^^^ = 9 x 10~^^ which might become 



possible result of the Pathfinder mission 1 11 511. 

Let us now consider some generalizations. In ^ f{fV^'^) has the form 
Consider now the minimal modification of the form 



f{rlV') = {rlV'Y-^. (16) 

According to (P) this modification of the graviton propagator gives the fractional 
change of the gravitational potential analogous to (3): 

e ~ r;i-%-i/V3/2+^. (17) 

The corresponding anomalous lunar perihelion precession is 

50 ^ r;i-V;i/V=^/2+7, (18) 



^We thank these authors for sharing their prehminary results with us. 
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Observations of accuracy can therefore test gravity theories with 



/ r \ ^(^+^' 

One can speculate that in the absence of additional scales the gravity theories 
that produce self-acceleration, without vacuum energy should have Vc ~few Gpc. 
Then the lunar precession accuracy of ~ 10^^^ will tests the 7 = theory [0. 
The dependence of the right-hand side of (|^) on 7 is very strong, and cosmologically 
interesting theories with 7 < are ruled out by current observations, while theories 
with 7 > are not testable by the solar system observations. 



4 vDVZ Discontinuity and its Absence 

Some time ago van Dam and Veltman, and Zakharov suggested that the solar 
system observations rule out the possibility of non-zero graviton mass, no matter 
how small. Their conclusion was based on a linearized theory of massive graviton 
with the following action 



2' 



S = Sei+ I d'x ( - {hff) + '-h.^T^'^ 1 (20) 



The first term on the r.h.s. is the standard Einstein action expanded to a quadratic 
order in the metric fluctuations about the flat space g^^ = rj^^, + h^^. The mass 



term in (^), has the Pauli-Fierz form |2^, which is the only possible ghost-free 
combination quadratic in h^^. 

Discontinuity is due to the fact that massive graviton contains five degrees of 
freedom (five polarizations) as opposed to two polarizations in the massless case. 
Three out of the five polarizations couple to the conserved energy momentum source, 
leading to an additional scalar attraction at distances r -C as compared to 
the massless case. Since the additional degree of freedom couples differently to 
relativistic and non-relativistic sources, the effect is not merely reducible to the 
rescaling of the Newton's constant Gn-, and is observable at the level of one-graviton 
exchange. Thus, at this level the theory is discontinuous in the limit rUg 0. 

The amplitude of the lowest tree-level exchange by a single massless graviton 
between two sources with energy-momentum tensors T^^, and T^^ is (the tilde sign 
denotes the quantities which are Fourier transformed to momentum space): 

In the massive case this amplitude takes the form: 



•^massless ^2 \'^^^'^ r) j . (21) 



_ S-kGn _1 fp\f,^u .22) 

•^massive o , n I ui' o H'-'^ d I " V / 

g^' + mjV 3 ^ J 
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The additional scalar attraction for massive case, is reflected in the difference in 
the tensor structure. This difference can not be eliminated by simple redefinition of 
parameters and is finite, for arbitrarily small graviton mass. 

This result goes under the name of vDVZ discontinuity, and if true, would be an 
extremely powerful result, as it would rule out not only the possibility of massive 
gravity, but much wider classes of theories that modify gravity in far-infrared. 

However, the story is not so straightforward as we shall now discuss. Shortly 
after vDVZ observation, Vainshtein |^ suggested that vDVZ discontinuity was an 



artifact of the linearized approximation, and would be absent in fully non-linear 
theory. That is, vDVZ result was obtained at one-graviton exchange level, that is in 
the first order in expansion. This corresponds to solving the linearized equation 

+ (/.^^ - v.uh'^J = (23) 

However, Vainshtein noted that perturbative expansion m Gn breaks down in the 
zero TUg limit. 

However, he also showed that if one takes into account non-linearities, then for 
small TUg the perturbative expansion in powers of nig can be organized, in which 
case discontinuity can disappear. Since no fully non-linear generally-covariant theory 
for massive graviton was known, Vainshtein used a theory obtained by non-linear 
completion of only the first, Einstein, term in r.h.s. of (|2^) : 

G^, + ml{h^, - 7]^,h2) = SttGjv T^u (24) 

where Gf^u is the Einstein tensor. In which case solving for goo = e'^^'^^ and grr = 
components of the spherically-symmetric metric outside massive body, he found: 

z/(r) = — — -I- O {rn^g\Jrgr^^ , A(r) = — + O (^^g\JV~r' 

(25) 

Note that in the same parametrization the standard Schwarzschild solution of the 
massless theory takes the following form: 



(26) 

Here Tg = 2GnM is the gravitational radius of the source of mass M. 

This expression is fully continuous in mg and reproduces Einsteinian results in 
the zero- mass limit. Thus vDVZ has indeed disappeared at the non-linear level. 
Notice the two important facts. First the sub-leading corrections to the metric are 
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non- analytic in Vg. Secondly, the deviations from the standard Einstein gravity 
become important at distance 

r. = ^^^^, (27) 
rrig 

which is parametricaly shorter than the Compton wave-length of graviton. This is 
not surprising, if we recall that vDVZ was cured by non-linearities that are most 
important in the neighborhood of the heavy sources (large r^). As a result the heavier 
is the source, the larger is the critical distance at which linear approximation takes 
over. 

The above results however, may not sound completely satisfactory, since the 
mass term in (p^ was still kept at the linearized level, and as a result theory was 
not fully generally-covariant. As we shall see, the above results nevertheless persists 
in theories, which are fully non-linear and generally-covariant. There too the dis- 
continuity is absent and corrections can penetrate at distances much shorter than 
the scale at which linearized gravity gets modified. The latter fact gives possibil- 
ity of experimentally testing such theories through precision measurements at the 
relatively short distances. 

An example of generally-covariant theory that modifies gravity in far infrared, 
and which exhibits vDVZ at the linearized level is given by the following action 

/\/f2 p I i\/r2 „ — 

S = ^ d^xdy 7^ + ^ d^x^/\g\ R{x) . (28) 

J 2 J * 

Where is 5D metric tensor, A, B are five- dimensional indexes, and TZ is the 
five-dimensional Ricci scalar, g^^ denotes the induced metric on the brane which we 
take as 

9A^) =gtK^.y = ^), /i,z/ = 0,l,2,3, (29) 

neglecting the brane fluctuations. 

We assume that our observable 4D world is confined to a brane which is located 
at the point ?/ = in extra fifth dimension. That is the energy-momentum tensor 
of 4D matter has the form T^^{x) S{y). 

Although, the underlying theory is high-dimensional, from the point of view of 
a 4D observer localized on the brane, Newtonian gravitational potential is just the 
usual 1/r gravity, which gets modified to only at very large distances r 3> Vc- 

The gravitational potential between the two bodies located on the brane can be 
read-off from the form of the Greens function (it is convenient to work in momentum 
space in the four world-volume directions and in position space with respect to the 
transverse coordinate y). 

Neglecting the tensorial structure of the propagator the scalar part of the Green 
function has the following form 0: 

G{q,y = 0) = ^ (30) 
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For the static gravitational potential \E'(r), one gets at short distances, i.e., when 



TT 



+ 



-1 + 7 - In 



(31) 



Here 7 ~ 0.577 is the Euler constant. The leading term in this expression has the 
familiar 1/r scaling of the four-dimensional Newton law. For r ^ one finds: 



•qj(r) 



1 



levr^Mf r 



1 / 1 

+ O - 



(32) 



The long distance potential scales as 1/r^ in accordance with the 5D Newton law. 

The above model (^) exhibits the vDVZ discontinuity in the one-graviton tree- 
level exchange. This can be seen directly from the /zz/ components of linearized 
Einstein equation for bulk metric fluctuations about the flat 5D metric, which after 
gauge fixing can be brought into the following form B 



6{y) + 5{y)d^d,hl .(33) 



The first term on r.h.s of this equation has a structure which is identical to that 
of a massive 4D graviton. The second term, d^dy, vanishes vanishes whenever it is 
contracted with the conserved energy-momentum tensor, and thus plays no role at 
one graviton exchange level. As a result, the amplitude of interaction of two test 
sources takes the form: 



oc 



1 n^pt-n^'f 
3 -'m 



(34) 



where q = \fcp'. We see that the tensor structure is the same as in the case of the 
massive 4D theory, see Eq. (p2D, which signals vDVZ discontinuity. 

The discontinuity however is cured at the nonlinear level. This was demonstrated 



explicitly both for cosmological solution ||Tl 
Schwarzchild 



as well as for cosmic strings |12| and 



131, 161. We shall concentrate on the latter case. The solution of 13 



has the form: 



i/(r) 



r 



A(r) 



(35) 



This expression exhibits some interesting behavior. First, we see that the two fea- 
tures observed by Vainshtein persist here, but there are some peculiarities. The 
correction is non-analytic in r^, and moreover it becomes important at distances 
much closer than Tc. This is in accordance with Vainshtein case. However, there 
are important differences which are crucial for observations. Notice that, one would 
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naively expect that if Vainshtein's analysis is correct, then the above expression 
should be obtained from Vaishtein case by substitution mg — >• l/vc- The reason is 
that, although we are not dealing with massive gravity in the strict sense, never- 
theless, Tc is the scale at which 1/r-law of linearized approximation breaks down. 
It is true that in our case it is replaced by rather than by exponentially sup- 
pressed Youkawa potential, but this difference naively seems inessential. So naively 
one would conclude that Tc should play the role of mg in controlling the strength 



of correction. However, this is not the case. Instead, the expression of ref [0 is 
obtained from Vainshtein by substitution: 



1 



m^^ — (36) 

We shall now try to understand why this is the case, but let us first explain 
non-analyticity of the correction in Vg. The above theory contains two parameters 
Mpi and r^. vDVZ happens at the linear level due to the term that is proportional 
to I/tc, and is cured by the non-linear corrections proportional to Mp^. These 
non-linear corrections make sure that metric gets standard Einstein form near the 
sources. The deviation from the standard metric is due to the fact that linearized 
approximation becomes good again and 1/rc-terms take over. Thus, corrections to 
the metric are due to the fact that one-particle exchange dominates at large scales. 
Then the leading corrections should be proportional to 1/ Vc- Close to the source the 
metric can be expanded in powers of I/tc, e.g. 

T 1 

5'oo('^) ^c, ^g) = 1 + — -\ ^a^^-a _|_ j-^igher order terms (37) 

r Tc 

where a is some positive power. Now, we know that for the fixed r, for weak 
sources, the linearized approximation is valid, and thus the effect of extra graviton 
polarization must be reintroduced. Thus, for fixed r, the second term on r.h.s. of 
(0) must dominate in the limit — 0. This can only happen if a < 1, which 
explains non-analyticity of the correction. To summarize shortly, given the fact 
that difference in the tensor structure of the metric is introduced by dominance of 
one-particle exchanges, which are proportional to ^, and dominate for small r^, the 
corrections to Schwarzchild cannot be analytic in r^. 

To understand which quantity plays the role of m^, the following observation is 
useful. If we are interested in making 4D metric on the brane continuously approach 
Schwarzchild in the limit 1/rc — 0, it is enough to have non-linear action on the 
brane only, and keep high- dimensional part of the action linear. In other words 
for curing vDVZ the non-linear interactions in the bulk are unimportant. That is, 
for obtaining Schwarzchild solution in the leading order on the brane the following 
equation is enough 

Mil {^0^. + Ky)G,u ) = 5{y)T,, (38) 
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Where, Qj^^ is the five-dimensional Einstein tensor, hnearized on a fiat background. 
This system describes the five dimensional graviton h^^{x,y) that freely propagates 
in the bulk, and has non-linear self-couplings only on the brane. These brane- 
localized self couplings are most important for determining the metric on the brane 
near the source localized on the brane. Diagrammatically, this fact can be un- 
derstood as follows. The nonlinear corrections that cure vDVZ discontinuity cor- 
respond to all the tree-level diagrams with virtual graviton lines that end on the 
sources localized on the brane. The self-interaction vertexes of virtual gravitons can 
be located both on the brane or in the bulk. The bulk non-linearities come from 
the 1/rc-suppressed bulk curvature term and are subleading. So the diagrams for 
which graviton self-interaction vertices are not located on the brane are sub-leading 
in 1/rc-expansion. Theory with linearized bulk action gives rise only to the di- 
agrams in which all graviton vertices are located on the brane. Since these diagrams 
are dominant, the eq ( |55D is suffices to cure vDVZ discontinuity. Now let us compare 
the five-dimensional theory defined in by eq. (|38|) to an effective four- dimensional 
theory of four-dimensional graviton h^j,{x) defined by the following equation 

Mil [g^, + ^v^(V - ri,uK)^ = T^u (39) 

These can be treated as two independent theories, but some sub-class of amplitudes 
in the two cases are equal. In fact the interaction amplitudes between brane-localized 
sources with no graviton emission in 5D theory are equal to the ones of 4D one. 
Notice that the two theories are designed in such a way, that the Greens function of 
the four-dimensional theory, is equal to the greens function of the five- dimensional 
one evaluated at the point y = 0: 

{hf,y{x)Kp{x')) = {h^^{x, y = 0)hap{x', y = 0)) (40) 

Also all the nonlinear interactions in two theories are the same. Given these facts, 
it is obvious that any interaction amplitude between the brane-localized sources 
Tfj,u{xi), Tf^i,{x2), ... T^jy(x„), (in which no gravitons are emitted in the final state) in 
five dimensional theory (|38|) will be equal to a similar amplitude in four-dimensional 
theory (0). For instance consider the lowest three- level interaction among the three 
sources with a single intermediate three-graviton vertex V{x)6{y): 

j d^xidyid'x2dy2d''x3dy3d'xdyT''%xi)6{yi) T"^ (0:2)5(2/2) T^'ix^) S{ys) 
y/.V',a'/3',7V'(a,) ( ^1) hf,'u'{x, y) ) ( V(x2, 2/2 ) K'^'{x, y) ) 

{h^p{x-i, ys) /iyp/(x, y) ) = 

J d^xid^X2d^X3d^xT^"'{xi) r"^(x2) T''P{x3) V^M''^'.°'/3''7'p'(a;) ^ hp^ixi) hp,^,{x) ) 

{ hap{x2) Kip>{x) ) { hypixs) h^>p>{x) ) (41) 

The last expression represents a similar amplitude evaluated in theory of eq. (^). 
Generalization of the above relation to more complicated diagrams with arbitrary 
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number of internal vertexes and external sources is trivial. As a result any such 
amplitude in higli-dimentional theory eq. (|3^) has an equal counterpart in four- 
dimensional one. Thus, if for example, we are interested in the metric created 
by the brane-localized source on the brane in theory defined by eq. (38), we can 
instead solve eq. (^) and get the correct answer. But for solving perturbatively for 
Schwarzschild metric, the equation (^) in the leading order is equivalent to that of 
Vainshtein (|2^) in which is substituted by l/rrg. 

Now we can generalize this result and give a simple qualitative prescription that 
gives a possibility to estimate the sub-leading correction to Schwarzschild in other 
theories of interest, in which Einstein gravity is modified in far infrared by adding 
some operators that at the linearized level behave as 



G,. + -/(r,^V^)(V - V,. K) = (42) 

^' c 

Then subleading correction to the metric can be estimated from (pSf ) by substitution 

(43) 




where action of the operator should be understood in terms of the eigenvalues of V. 
For instance: 

firlV') - = / rf'p/(ry) V (44) 
r J 

It should be noted that we have neglected corrections of the same order in l/vc 
coming from the nonlinear completion of the mass term. Thus our results is just an 
order of magnitude estimate. In principle there could be theories in which the addi- 
tional terms exactly cancel the corrections to Schwarzschild we computed although 
this is not the case in the theories of 0,0]. 



5 Applicability 



We shall now briefly formulate the applicability range of our analysis. Consider a 
theory which, in some gauge, at the linearized level reduces to eq. (^. The criteria 
for the presence of discontinuity (/5 7^ 1) is the existence of the following spectral 
representation: 

1 r°° dsp{s) 



s + 



(45) 



with a semi-positive definite spectral function p(s), which is requirement for the 
absence of unphysical negative norm states. In such a case the equation (|^) is 
derivable from the action 



S 



f 

Jo 



ds 



S 



is) 
El 



d X 



+ 



2M, 



PI 



(46) 
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where /i^,^ = / ds^p{s)h^^^. This action describes the hnearized theory of continuum 
of free massive gravitons with masses = s. They all couple to the same conserved 

source but with mass-dependent coupling -2^~- Therefore the propagator of this 
theory is simply an integral over the continuum of the massive propagators taken 
with the weights p{s) 

A{q) = -MgJ dsp{s)^ ^\ ' (47) 

JO s + 

where in our parametrization c(0) = 3/2 and c(s 7^ 0) = 1. 

A useful illustrative example, in which vDVZ can be seen from the above rea- 
soning is the model of 0. There the representation has a simple physical 
meaning It is just an expansion into the continuum of massive Kaluza-Klein 
states. In that case the spectral function is given by p{s) = ^^^^^2^^ . In this 
language the presence of vDVZ discontinuity in is trivially understood. Since 
each Kaluza-Klein mode is a massive spin-2 particle, each of them exhibits vDVZ 
discontinuity by default. Thus exchange by tower of Kaluza-Klein states exhibits 
the same discontinuity as the individual exchanges do. 



6 Conclusion 

There are strong motivations coming from cosmology for modifying the standard 
Einstein gravity at large distances. A wide class of such gravity theories can be 
tested by astronomical observations of the solar system. 
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